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The major effort has been to understand the morphology of the ion cloud 
in the radial shaped charge barium injection experiments of Wescott et al 
(1980). Basically these experiments consisted of projecting barium 
perpendicular to the earth's magnetic field by a radially symmetric explosive 
shaped charge. The experiments were conducted from rocket-borne payloads at 
an altitude where interaction with the ambient neutral atmosphere is 
negligible. ■ 

The bdoic feature we have attempted to understand is the shape of the ion 
cloud that remains after the explosive products and neutral barium clears 
away. Panel D of Figure 1 shows the ion cloud which has the configuration of 
a rimless wagon wheel. The major features are the 2.5 km radius "black hole" 
In the center of the cloud, the surrounding ring of barium ion and the spokes 
of barium ionization radiating away from the center. The cloud shows no 
evolution after it emerges from the neutral debris, so it must have formed 
within 5 seconds of the event. 

Panels A-C of Figure 1 show the expanding neutral barium. The most 
notable feature is the inhomogeneous or blobby structure of the expanding 
neutrals. The spokes in the ion cloud can then be explained as the trail of 
ionization left behind by a neutral blob as it photoionlzes. The ambient 
magnetic field is effective in halting the outward expansion of the ions as 
they are created. 

The existence of the black hole, however, requires the outward transport 
of Ionization within 2.5 km of the explosion, and the nature of the transport 


mechanism Is the major problem addressed In this work. To this end a 
numerical model was set up to calculate the motion of Ions and electrons 
subject to the electrostatic and Lorentz forces, and subject to being 
col 1 1 si onal ly driven by the expanding neutrals. The velocity of the neutrals 
was calculated based on the motion of neutrals deduced from television imaging 
of the event. The electric potential was determined from the quasi neutrality 
condition derived from requiring the divergence of the current to be zero. 

This gives an equation for the potential of the type 

V«EpV4>^7 2: H XV<t> = S (1) 

where z and e u are the field line Integrated Pedersen and Hall conduct- 

P H 

ivities respectively, and S is a source term proportional to the velocity of 
the barium neutral s. 

The investigation has been frustrated by our inability to find a 
convergent iterative scheme to obtain solution to (1). The large Hall term 
destabilized the iterative schemes we have tried. 

As an alternative we have modified the calculation by dividing the space 
into pie-slice segments of angular width a e where a© 3 2 n/n, n an 
integer. Each segnent was divided into two subsegments of angular width 
A©j + a ©2 where A 0 ^+ a ©2 3 A0- To represent the azimuthal variation in the 
density of neutrals, one neutral gas density was assigned to the segment a 0 j 
and another to the a ©2 segment. The pattern repeats n times. The grid lines 
are on the lines separating the aQj and A0 2 sectors. The potential is 
determined by the condition that the current remain divergence free. 

Electrons can flow from one sector to another, but the ions are assumed to 
have no azlnuthal drift. 

The result from this simplified model is that the polarization electric 
field allows the uninpeded transport of plasma across the magnetic field at 
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the velocity of the neutrals. Figure 2 shows the computed Ionization density 
at 1 sec after the event for the case where the neutral density In the 
A0^ and a ©2 segments was the same. This was at a time when all Ion motion 
had ceased. The large peak In density near the origin can be seen. This is 
Inconsistent with the existence of the black hole in Frame D of Figure 1. 
Figure 3 shows the results when the density In the two sectors differed by a 
factor of 10. The density minimum at the origin is consistent with the black 
hole. 

The calculations show that the azimuthal variation in the density of the 
neutrals is critical to the formation of the black hole. The limitation of 
the model is that it does not allow azimuthal motion of ions, which we believe 
is important to the formation of the ring of barium ion around the black 
hole. The results of this work were presented to the Fall 1984 meeting of the 
American Geophysical Union in San Francisco (Hollerbach et al , 1984). 

A more definltve simulation of the radia’ shaped charge experiment will 
require the solution to (1). We are proceeding along two tracks. One is to 
calculate the number density by transporting the electrons rather than the 
ions. This should produce density gradients more nearly parallel to the 
potential gradient and thus reduce the magnitude of the destabilizing form in 
(1), ve h x v<f> . If this does not result in convergent iterations, we shall 
use sparse matrix methods to obtain a direct solution to (1). 
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Figure Captions 

Figure 1. A sequence of TV pictures showing the evolution of the King Crab 
radial barium injection experiment (Wescott et al . , 1980). 

Panel A and B show barium neutrals, while in C a ring of barium 
neutrals can be seen expanding away from a faint barium ion 
cloud. Panel D shows the barium ion cloud through a filter. 

Figure 2. Calculated radial distribution of barium ion density at one 
second computed for the neutral velocity distribution in the 
Bubble Machine experiment, assuming a radially symmetric 
expansion. Hie motion of the ions is stopped by the electrons 
which remain rigidly tied to magnetic field lines. 


Figure 3. Similar to Figure 2, except that the time is at 3 seconds and 
that a 10 to 1 azinuthal variation in neutral density is 
assuued. The barium is moving at essentially the neutral 
expansion speed. 
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